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Abstract 

We show that a discrete exchange symmetry can give rise to realistic dark matter 
candidates in models with warped extra dimensions. We show how to realize our 
construction in a variety of models with warped extra dimensions and study in 
detail a realistic model of Gauge-Higgs Unification/composite Higgs in which the 
observed amount of dark matter is naturally reproduced. In this model, a realistic 
pattern of electroweak symmetry breaking typically occurs in a region of param- 
eter space in which the fit to the electroweak precision observables improves, the 
Higgs is heavier than the experimental bound and new light quark resonances are 
predicted. We also quantify the fine-tuning of such scenarios, and discuss in which 
sense Gauge-Higgs Unification models result in a natural theory of electroweak 
symmetry breaking. 



1 Introduction 



Models with warped extra dimensions [T] have arisen in the last few years as strong candidates 
for a natural theory of electroweak symmetry breaking (EWSB). The original solution to the 
hierarchy problem has been supplemented by the addition of a natural flavor structure [2] , 
a custodial symmetry to protect the T parameter [3] and the Zb^bi coupling [3], and the 
realization of the Higgs as the pseudo Goldstone Boson of a broken global symmetry [5]-[7]. 
These developments have produced calculable models that successfully address most of the 
mysteries related to the electroweak (EW) scale. 

On the other hand, dark matter (DM), that is often considered one of the strongest - 
albeit indirect- hints of physics beyond the Standard Model (SM), has so far lacked a generic 
implementation in models with warped extra dimensions. The main reason is the inherent 
asymmetry in warped backgrounds, that do not posses a natural KK parity as the one present 
in Universal Extra Dimensions (UED) [8]. 

In this article we explore a generic procedure to introduce an exact discrete exchange 
symmetry that results in new stable states, without introducing new parameters. This is 
done via a doubling of part of the field content. The exchange symmetry we advocate has 
been first introduced in [S] (where it was dubbed "mirror symmetry") and [TU] to alleviate 
the fine-tuning and to get a viable DM candidate in Gauge-Higgs Unification (GHU) models 
in flat space. As already anticipated in [9], it can be extended straightforwardly to warped 
models. Given a bulk field, <j>, satisfying certain boundary conditions (b.c), the procedure 
consists of replacing cf) by a pair of fields, fa and fa, and imposing the symmetry fa <-> fa. []] 
The even linear combination <ft+ = {fa + fa)/\^2 is identified with the original field (in 
particular, it inherits the b.c. obeyed by 0, as well as its couplings). The couplings of the 
orthogonal combination, = {fa — fa)/\/2 are determined by those of Under the 
exchange symmetry one has <f>± — ► ±(j)±, so that one can assign a multiplicative charge +1 
to (j) + and —1 to Provided the discrete exchange symmetry is an exact symmetry at 
the quantum level, the lightest Kaluza-Klein (KK) resonance among all Z2-odd states in the 
model is absolutely stable, and will be referred to as the LOP (Lightest Odd Particle). We 
argue that the above symmetry is indeed exact, by showing that possible 5D Chern-Simons 
(CS) terms, in general needed to restore gauge invariance in 5D theories, do not violate it |12j . 

The choice of which fields to double must be guided by phenomenological considerations. 
For example, DM direct and indirect searches impose stringent constraints on the possible 
couplings of the DM candidate to SM fields. It is therefore natural to look for charge and 
color neutral fields that can lead to viable DM candidates. In fact, the models we will 
consider always contain U{\) factors, and it will be natural for the DM candidate to be a 

A similar discrete symmetry has recently been used in [11] to get a DM candidate in the context of little 
Higgs theories. 
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U(l) massive gauge field, X_. This is similar to the 5D UED case and the GHU model of [10] 
in flat space, in which the DM can be identified with the first KK mode of the hypercharge 
gauge field. If the LOP were the only Z2-odd particle, it would couple to SM fields only 
via non-renormalizable interactions. For the typical scales involved, its annihilation rate 
would then be extremely small, and the resulting thermal relic density unacceptably large. 
Hence, it is necessary in general to apply the previously described "doubling" construction 
to additional fields, making sure that X- remains the lightest among the Z2-odd particles. 
We can choose b.c. for the LOP so that it has a mass about one order of magnitude smaller 
than the mass of the first Z2-even gauge resonances. Thus, our DM candidate has typically 
sub-TeV masses, which, as we will see, can range from ~ 100 GeV in Higgsless models up to 
~ 700 GeV in Randall-Sundrum (RS) models. Together with its mass, the coupling of the 
LOP to SM matter is the other crucial parameter governing its relic density. In all the cases 
we will consider, the U(l) symmetry is related to the SM U(1)y symmetry, with a coupling 
constant of electroweak size, and the Z2-odd fields will consist of X_ and a subset of fermion 
fields, typically associated with the top or bottom quarks. It is clear that our DM candidate 
is a weakly interacting massive particle (WIMP) with approximately the correct mass and 
couplings to give rise to the observed DM relic density in the universe. 

The construction outlined above turns out to be particularly natural in the specific class 
of GHU scenarios in warped spacej^] In this case, the additional set of fields introduced by 
our procedure is not only minimal, but is also singled out by the role the doubled sector plays 
in generating the EW scale dynamically. Although it is not easy to use the discrete sym- 
metry to protect EW observables from large tree-level corrections, as in R-parity preserving 
supersymmetric scenarios, Little Higgs theories with T-parity, or UED with KK parity, we 
will see that the physics of EWSB in the GHU model with the discrete symmetry naturally 
leads to a region in parameter space where the constraints due to precision measurements are 
relaxed. Furthermore, this same region of parameter space, in which coannihilation effects 
can be relevant, predicts a dark matter abundance in accord with observation. Thus, the 
DM sector is tightly connected to the physics of EWSB, and plays an indirect role in leading 
to agreement with precision constraints J^] Also, the model is rather predictive, with several 
fermionic resonances nearly degenerate with the DM particle that should lead to an exciting 
collider phenomenology. 

We also consider some non-perturbative effects -the formation of bound states and the 
effect of QCD Coulomb-like forces- which might invalidate the usual perturbative computa- 
tion of the relic density. We argue that the problem of the formation of bound states does 

2 It is known that such models can also be seen, thanks to the AdS/CFT dictionary [13], as strongly coupled 
4D composite Higgs models. However, since calculability requires the 5D picture, we will mostly use the 5D 
language, adopting only occasionally the 4D dual language. 

3 See [T3] for a model in which EWSB and DM are related. In this reference, however, the hierarchy problem 
is not addressed. 
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not occur and that the effect of Coulomb-like forces on the perturbative cross-sections is 
negligible. 

We also discuss in some detail the degree offine-tuning involved in the EWSB pattern, and 
compare to supersymmetric extensions of the standard model. It turns out that in the GHU 
framework (with or without DM) the fine-tuning is somewhat worse than expected by naive 
considerations. We point out, however, that here the fine-tuning seems to be associated with 
accommodating a top mass of order the EW scale, rather than with an intrinsic tension in 
the Higgs sector. When restricted to the region of parameter space with a fixed top mass, the 
low-energy properties of the model turn out to be essentially insensitive to the microscopic 
parameters of the model, and therefore to the detailed properties of the new physics, which 
is a very interesting feature. 

Although we find the application of our construction to the GHU framework particularly 
appealing, we stress that it is of more general applicability. We also briefly discuss the 
implementation of the discrete exchange symmetry in two other models with warped extra 
dimensions: the simplest RS model with the SM fermions and gauge bosons in the bulk, and 
a Higgsless model. Hence, our construction leads to viable DM candidates in a variety of 
scenarios, without the introduction of new parameters. 

The organization of the paper is as follows. In section[2]we review the generic properties of 
the exchange Z2 symmetry that gives rise to a dark matter candidate. In section [3] we discuss 
in detail its implementation in a model of GHU. In particular, we study the interplay between 
EWSB, the EW constraints and the DM relic abundance, the (ir)relevance of the above 
mentioned non-perturbative effects, and discuss the fine-tuning. In section [4] we describe the 
implementation of the discrete exchange symmetry in other models. Section [5] is devoted to 
particular details of DM collider phenomenology and DM direct detection in our construction. 
We comment on the issue of anomalies in section [6j and conclude in section [7j We relegate 
some technical details to the Appendices. 

2 Exchange Z 2 Symmetry and Stable Particles 

The simplest way to get a stable massive particle X- beyond the SM spectrum is obtained 
by introducing a discrete Z2 symmetry under which all SM particles are even, with X_ 
the lightest Z2-odd particle. For example, viable DM candidates in promising models of 
new physics, such as SUSY or little Higgs models, are stable due to Z2 symmetries (R- 
and T- parity |15| . respectively). In the context of extra-dimensional theories, a geometric 
Z2 symmetry called KK-parity leads to viable DM particles in 5D [16J and 6D [17J UED 
scenarios. Most extra-dimensional models which aim at stabilizing the electroweak scale are 
based on warped compactifications of the Randall-Sundrum type [T] , with bulk gauge and 
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fermion fields [2] , where KK-parity cannot be trivially imposed^] We are therefore motivated 
to look for other discrete symmetries that could lead to DM candidates in such warped 
scenarios. 

From a model-building point of view, it is always possible to impose a discrete symmetry 
ad hoc for the only purpose of getting a stable particle, possibly with the correct properties 
to account for the observed DM density. Such constructions become far more interesting if 
the DM sector is tightly connected with other sectors of a given theory, in such a way that 
it leads to additional testable predictions and/or no new parameters are introduced. The Z2 
exchange symmetry we consider belongs to this class, since no new parameters are added. 
Furthermore, in the particular class of GHU models, which will be considered in detail in 
the next section, it can also improve the EWSB pattern and lead to better agreement with 
electroweak constraints. 

Consider a model with warped extra dimensions and a bulk gauge symmetry that includes 
a U{\) factor, Q x U{\)x- Our construction introduces a discrete symmetry that acts on the 
U(l)x factor and leads to a stable spin-1 particle. The group Q does not play any role in our 
construction (it will be neutral under the discrete symmetry) and we will not discuss it any 
further. The first step is to double the U{l)x field, U{l)x — * U(l)x 1 x U(l)x 2 - The original 
U(l)x gauge boson is identified with the symmetric combination, X + = (X\ + X2j/\/2, and 
the gauge couplings chosen as gx 1 = 9x 2 = V%9X > so that X + couples with strength gx • The 
antisymmetric combination, X_ = (X\ — X2)/\2 is odd under the Z2 exchange symmetry 
under which X\ *-* X2. 

As remarked in the introduction, doubling only a U(l) symmetry is not enough to obtain 
a realistic DM candidate. In order to improve this, we take a subsetj^] of the fermionic 
fields in the model ip, with U(\)x charge Qx, and double them into mirror pairs, ip — ► ip\2, 
assigning them U(l)x 1 x U(l)x 2 charges (Qx, 0) and (0, Qx), respectively. Finally, undoubled 
fermions, cp, with charge Q under the original U(l)x group, are assigned charges (\Q, \Q), 
which implies that they couple only to the symmetric combination X + , with charge Q. Since 
ip\ couples only to X\, while 1P2 couples only to X2, the relevant 5D Lagrangian density is 
C = £1 + £2 with 



e- 1 b = ~F? + ilH 



h M {V M + iQxgXiXiM) -miil>i, i = 1, 2 , (2.1 



where Fi is the field strength of the gauge field Xi with 5D coupling constant gx t , T^M is the 
gravitationally-covariant derivative and e is the determinant of the funfbein associated to the 

4 Ref. |18j considered a discrete symmetry to forbid low scale baryon number violation in Grand Unified 
Models, that also gave rise to a stable DM candidate. See also [19] for a recent attempt to include KK parity 
in warped space, and [20] for some proposals of DM candidates in composite Higgs models. 

5 Doubling the whole spectrum is unnatural, since fermions with twisted boundary conditions can have 
exponentially small masses [211 118] . giving rise to unacceptable LOP's. In fact, as we will see, requiring X- 
to be the LOP constrains |c| < 1/2 for the doubled fermions. 
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5D metric 



ds 2 = e- 2ky 7]^dx^dx u -dy 2 =(—y (r ttiv dx^dx u - dz 2 ) , (2.2) 
where, as usual, \x runs over the 4D directions, < y < L, z vv < z < z m , z = e ky /k. In 
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Eq. (2.1), we have not explicitly written terms containing the gauge bosons of the group Q 



or the undoubled fermions. The exchange symmetry constraints m\ = m<i = m = ck. Notice 



that for simplicity we omit a possible mixing term of the form F1F2 in Eq. (2.1). Similarly, 
depending on the b.c. chosen for the Z2-odd fields, possible boundary terms can appear. 
We assume here that all these operators can be neglected, so that no new parameters are 
introduced. In terms of ± fields, we have 

e- x L = ~(F 2 + F 2 ) + Mi^Vm - m)%b + + ^_(i 7 M ^M - m)^- 

- gxQxl M $+X +M ^+ + 4>-X +M ip- + $+X- M *l>- + ^-X_ M ^+) ■ (2.3) 

Let us now discuss the boundary conditions for the different fields. As explained above, X+ 
inherits the b.c. of the original U(l)x gauge boson (that can in general involve in a non- 
trivial way the neutral gauge bosons in the group Q). The b.c. for the odd combination are 
taken to be 

X_ li = ±={X Xll -X 2ll ) ~ (+,-) , (2.4) 

where +/— denote Neumann/Dirichlet b.c, respectively, with the first/second entry in paren- 
thesis referring to the UV/IR boundary (X_ 5 satisfies opposite b.c). These b.c. allow for UV 
brane localized kinetic terms characterized by a dimensionful coefficient r^ v (in the notation 
of Ref. |22]). The lightest X_ resonance (the LOP) has mass of order 



mx -"V^^) MlR "v^ Mm ' (2 ' 5) 

where fi m = ke~ kL = l/z m and the second equality holds whenever the localized term is 
small. This mass is parametrically smaller than the KK scale For example, for values of 
kL that solve the hierarchy problem, mx_ is about a factor of 10 below the mass of other 
gauge resonances, which are of order mx + — 2.5 /i IR . Note that (— , +) b.c. for X_ would 
instead give a larger mass, of order mx + - It would then be hard to identify X_ as the LOP 
and get the correct relic density with such a choice of b.c, without the introduction of large 
brane kinetic terms. For these reasons, we do not consider this possibility. 

Regarding the doubled fermions, satisfies the boundary conditions of the original 
fermion. The odd combination, ^_ should obey (+, — ) or (— ,+) b.c. (for one chirality, 
opposite for the other one) so that no fermion zero modes are introduced in the Z2-odd 
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sector of the theory]^] However, the unbroken U{1)- symmetry on the UV brane, that mixes 
the two mirror fermions ip±, requires that the b.c. of -0- be equal to that of tp+ on the UV 
brane. The mass of the first KK mode of ip- strongly depends on c. For |c| > 1/2, this 
lightest fermionic KK mode can be lighter than X_ . In the context of warped scenarios that 
solve the flavor puzzle by fermion localization, this usually means that we can only double 
the fermions associated with the top and bottom quarks, which have |c| < 1/2. 

These are the qualitative, model-independent features of our mechanism to endow existing 
models with a stable particle that can be a DM candidate. In the next section, we apply this 
construction to a very appealing model of EWSB in warped extra dimensions, studying in 
detail the pattern of EWSB, EW precision tests and the calculation of the DM relic density. 
In section [4] we will briefly show how this construction can be easily implemented in a variety 
of models with warped extra dimensions. 

3 A Relevant Case: a GHU/Composite Higgs Model 

Models of GHU (see e.g. [23J for a brief review and further references) rely on an extended 
symmetry Q broken to a subgroup TC. The KK gauge bosons, A^, associated with the coset 
directions Q /TC are all massive, but the A§ towers give rise to zero-modes. The latter are 4D 
scalars and, by an appropriate choice of the symmetry breaking pattern Q — > TC, can have 
the correct quantum numbers to be identified as the Higgs. The higher-dimensional gauge 
symmetry leaves a remnant shift symmetry for the Higgs [23] that ensures that the Higgs 
potential is finite to all orders. This UV insensitivity guarantees that the Higgs potential 
gets corrections of the order of the IR scale, fj> m , and not the cut-off of the theory, thus 
alleviating (though not solving) the little hierarchy problem. In this paper we will show 
how our construction can be easily implemented in a fully realistic model of GHU in warped 
extra dimensions. We will consider the minimal composite Higgs model with fermions in the 
fundamental representation of SO(5) (MCHM5) of Ref. [25], but our construction could be 
applied to other variations such as the models considered in [26J. The starting bulk gauge 
group is SO (5) x U(l)x, broken to 50(4) x U{l)x on the IR brane (the actual symmetry 
on the IR brane is assumed to be 0(4)) and to the SM on the UV brane. The quark 
sector is embedded in fundamental representations of SO (5) that decompose under 50(4) ~ 
SU(2)l x SU(2)r as a bidoublet plus a singlet, 5 = (2,2) © (1, 1). The relevant sector for 
EWSB is the third generation. See appendix [A] for a summary of the field content of the 
model. 

Let us now construct the simplest Z2 extension of the MCHM5 model. Following our 

6 In principle, it is also possible to introduce Z2-odd fermions with (++) boundary conditions, getting rid 
of the zero modes by coupling them through mass terms to localized chiral fermions. In the limit that these 
masses become large, a description by effective b.c. as the one we are using is appropriate. 
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prescription, we enlarge the gauge group by doubling the U(l)x factor to U(l)x 1 x U(l)x 2 
and identify the original gauge boson X with the even combination, X + = A^(Xi + X2). 
The odd combination has (+, — ) boundary conditions as in Eq. (2.4). As emphasized in 
the previous section, only fermions with mass parameter \c\ < 1/2 can be safely doubled, 
since otherwise they could easily give rise to charged/colored Z2-odd KK modes lighter than 
the Z2-odd KK gauge bosons. It is therefore natural to consider doubling the multiplets 
associated with the top quark, since accommodating the top mass requires |c| < 1/2 for the 
associated 5D fields. By looking at Eq. ( A. 5 ), it is clear that the only field that can be doubled 
without introducing new unwanted zero modes is £ u . We will see that this is sufficient to 
obtain the correct DM relic density. In summary, our quark sector is identical to that of 
the MCHM5 model, except that £ u is replaced by two copies, £ Ul) £ M2 , with U(l)x 1 x U(l)x 2 
charges (2/3, 0) and (0, 2/3), respectively, but otherwise identical in boundary conditions and 
bulk mass. The physical combinations are the Z2-even (which, for simplicity, we will call 
£ u in what follows) and the Z2-odd £~. The other quark and lepton fields are assigned the 
same charge under both abelian groups, (^Qx, \Qx), where Qx = 2/3 for Qx = — 1/3 
for £ 92 , and analogous assignments for the lepton sector. Therefore, the spectrum in our 
model contains a set of fields that are even under the Z2 symmetry and corresponds exactly 
to the spectrum in the original MCHM5 model, plus two 5D Z2-odd multiplets, X_ and 
The Z2-odd gauge boson has a first KK mode, which we will simply call X—, with a mass 
given by Eq. (2.5). The Z2-odd fermions have first KK modes with c u -dependent masses that 
are, neglecting EWSB effects and provided \c u \ < 1/2, always larger than the one of -X^-Jj 
As we take c u — > 1/2 the first KK mode of (2,2) u_ becomes degenerate with X_, whereas in 
the c u — > —1/2 limit it is (1, l) u_ that becomes degenerate with X_ (see appendix [A] for the 
notation) . 



3.1 Electroweak Symmetry Breaking and Precision Constraints 

The pattern of EWSB in the MCHM5 model and its main constraints have been studied 
in [251 [261 [27]. A realistic pattern of EWSB can be obtained for |c gi | < 0.4 and 0.35 < 
\cu\ i$ 0.45. Outside these two regions, it is difficult to obtain a reasonable value of the 
gauge boson, top and Higgs masses. Of particular interest is the dependence on c u . For 
fixed values of the other parameters, smaller values of \c u \ result in no EWSB, while larger 
values of \c u \ give the wrong EWSB pattern. The qualitative features of this dependence can 
be easily understood in the limit that the localized mixing masses vanish. In this simplified 
case, the two chiralities of the top quark arise from and obtaining a top Yukawa coupling 
of order one fixes c qi ~ 0.44. The contribution to the Higgs potential due to the £ 9l KK 

7 Note that odd-fermion masses are independent of the localized mixing masses, and are therefore entirely 
determined by c u , up to EWSB effects. After the mixing due to EWSB, these could become lower than mx_ 
if \c u \ is very close to 1/2. 
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tower destabilizes the origin, while the contributions from the gauge KK towers tend to 
align the vacuum along the EW symmetry preserving direction. It turns out that the top 
tower contribution is so large that it would drive the vacuum expectation value (VEV) to its 
maximum value, = 1, where 

f(h)\ 1 [Y . . 

8 h = sm{—) , h= g ]/kL^. (3.1) 

Here fh is the "Higgs decay constant" and g is the SM SU(2)l gauge coupling. The resulting 
EW symmetry breaking pattern is unacceptable since it leads to vanishing fermion masses, 
and highly non-linear couplings of the Higgs to the gauge bosons that are ruled out by 
EW precision measurements. However, fermions with "twisted" boundary conditions, such 
as £ u , give a contribution to the Higgs potential that tends to align the vacuum along the 
EW symmetry preserving direction. This contribution is controlled only by c u in the above 
simplified limit, and turns off when \c u \ > 1/2. If \c u \ <C 1/2, then this EW restoring 
contribution overwhelms the contribution due to £ 9l and results in Sh = 0. The upshot is 
that the desired EWSB VEV, < < 1, is obtained for \c u \ ~ 1/2, but not necessarily too 
close to 1/2. This also illustrates the crucial role that £ u plays for EWSB. When the mixing 
masses are turned on the details are more complicated, but the qualitative features remain 
the same. See appendix |B] for details on the computation of the Higgs potential. 

On the other hand, EW precision data typically prefer values of \c u \ close to 1/2. This is 
due to a sizable and negative one-loop contribution to the Peskin-Takeuchi [28] T parameter 
in most regions of parameter space [29], together with a non-negligible tree- level positive 
contribution to the S parameter. The tree-level contribution to S, when the light fermions 
are localized close to the UV brane, is given by [6] 



_, 67TS? 
'tree 



h - (3.2) 



g 2 kL nl 



IR 



where the second equality holds whenever v = (h) ~ 174 GeV <C fh- The negative contribu- 
tion to the T parameter arises at one- loop order from the lightest charge-2 /3 members of the 
SU(2)l x SU(2)r bidoublets, and is dominant away from the \c u \ ~ 1/2 region. Additionally, 
the pseudo-Goldstone nature of the Higgs leads to a further negative (positive) contribution 
to the T (S) parameter, due to the anomalous gauge couplings of the Higgs [30]. This can 
be described by an effective Higgs mass [31] 



m hteS = m h [ — , (3.3) 
where A is an effective cutoff scale of the order of the mass of the first SM gauge KK 
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resonances, ttikk, so that the corresponding shifts in S and T are 



1 . /"<eff N 



12-k \ ml I 
AT = 3 ln I — ^' eff 



Here w./i rcf is the reference Higgs mass used in the S-T fit to the EW data, and cw is the 
cosine of the Weinberg angle. 

The S-T analysis reveals that in the Z2 extended MCHM5 model, the regions with \c u \ ~ 
1/2 are preferred (see Fig. [TJ. For c u ~ —1/2, the reason can be traced back to the fact that 
the first KK excitation of the 5(9(4) singlet in £ u , that mixes with the top quark, becomes 
light and gives a positive contribution to the T parameter that compensates the negative 
contributions coming from the lightest bidoublet states and effective heavy Higgs, Eq. (3.3). 
As a result, one finds a sizable region in parameter space compatible with EWSB, where the 
EW precision measurements are relaxed. For c u ~ 1/2, it is again possible to satisfy the 
EW constraints, although the reason is somewhat more involved than for c u ~ —1/2. As 
emphasized in [29 , the charge-2/3 members of the bidoublets give rise to both positive (from 
the T| = +1/2, T% = -1/2 states) and negative (from the If = -1/2, T% = +1/2 states) 
contributions to T. In the custodially symmetric limit these two contributions would cancel 
exactly, but in most regions the fact that the T^ = —1/2 state is lighter than the T| = +1/2 
one, typically results in a net negative contribution to T (the signs are simply determined by 
the quantum numbers). However, when c u ~ 1/2, the lightest bidoublets arise mostly from 
£ u (and therefore have a nearly flat wavefunction component, that vanishes on the IR brane), 
and lead to a contribution to T via mass mixing with the bidoublet in £ 9l . It is then possible 
to suppress the coupling to the Higgs of the lighter Tj- = —1/2 state by ensuring that it lives 
mostly in the component with the nearly flat wavefunction, while the heavier = +1/2 
state has a larger component in the bidoublet of the £ 9l multiplet, that is localized near the 
IR brane. In this way, the positive contribution due to the T| = +1/2 state can dominate, 
and explains how the c u ~ 1/2 region can be compatible with the EW precision data with a 
relatively low scale 

In addition to the oblique corrections parametrized by S and T, non-oblique corrections 
associated with the third generation can also be relevant. These can be separated into flavor- 
preserving versus flavor- violating effects. A complete treatment of the latter would require 
the specification of the flavor structure of the model, which is beyond the scope of this paper. 
Most notable among the flavor-preserving effects are the corrections to the Zb^bi vertex, 
which has been measured at the few per mille level. Here we notice that the present model 



8 This important positive contribution is not present in the minimal models studied in Ref. [31] (it is, 
however, in models with a custodial protection of the Zb^bh coupling as the one we are considering). 
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enjoys the custodial protection pointed out in [4 , which reduces the tree-level corrections 
to this vertex to a level well below the experimental precision. However, as first pointed 
out in |29| [26] (see also [31 ) the one- loop contributions to this vertex can be significant 
and correlated with the one-loop contribution to the T parameter discussed above. Hence, 
we also include these effects in the fit to the EW precision observables, although, for the 
reasons explained above, we do not include flavor violating effects]^] In this analysis we have 
performed a fit to the Z-pole observables and the W mass [33] , and include the effects of four- 
fermion interactions that enter through the Fermi constant (these effects are subdominant). 
The upshot is that there is a bound of about fi m ~ 1.3 TeV in the c u > region, while the 
c u < region is somewhat more constrained, with a lower bound fi m ~ 1.7 TeV. 

Notice that although the Z2-odd fermions do not contribute to T (since the b.c. for 
£~ respect the custodial symmetry exactly and there is no mixing with custodial violating 
sectors of the theory due to the exact Z2 symmetry), they affect the minimization of the 
Higgs potential in an analogous way to £ M , discussed at the beginning of this section. In 
particular, the existence of these additional fields satisfying "twisted" b.c, and whose effects 
on the Higgs potential are controlled also by c u , means that the minimization of the potential 
favors values of \c u \ that are closer to 1/2 than if the Z2-odd fermions were absent. This 
is a welcome feature, since it goes in the direction preferred by the EW precision data, as 
discussed above. Furthermore, the opening of a well defined region with c u ~ 1/2, where the 
(tL,bL) doublet is more fundamental [see f^{z) in Eq. (C.2)] is also welcome since it might 
be argued, based on flavor considerations, that such a situation is more natural 

Thus, we find that the regions that lead to a correct EWSB pattern and good agreement 
with the EW precision data largely overlap in the Z2 extension of the MCHM5 model. We 
postpone a discussion of the details associated with Fig. [I] to the next section, after we have 
discussed the computation of the DM relic density in this scenario. For the moment, let us 
mention that we predict relatively light vector-like quarks (from both £ u and £ ~ ) with masses 
close to mx_ ~ fx m , a very distinctive signature of these models. It should also be 

emphasized that our model has the same number of parameters as the MCHM5 without DM. 



3.2 Calculation of the Relic Abundance 

As was discussed in the previous section, a realistic pattern of EWSB requires |c u | < 1/2. 
Thus, we have two separate regions that, a priori, could give rise to phenomenologically 
relevant scenarios. This can have important consequences for the computation of the DM 

9 Nevertheless, we have checked that a subset of the flavor-violating corrections, that are closely connected 
to parameters entering in the flavor-preserving effects [31], do not significantly alter the bounds. Specifically, 
these include certain loop- level contributions to the Zbs vertex that are bound by B —> X a l + l~ decays (these 
are likely to be more important in the c u < region), as well as tree-level contributions to the WtRbn vertex, 
that can be bound from the B — > X 3 y branching ratio [32 (these affect mostly the c u > region). 
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relic abundance since, as remarked above, some of the Z2-odd fermions become degenerate 
with the DM candidate when \c u \ = 1/2. In the case Cu ^ — 1/^5 ^ ^ the first KK mode of the 
50(4) singlet component of £~ that becomes close in mass to X- (the LOP), while the KK 
modes of the bidoublet are at least a factor of ten or so heavier. In the other relevant region 
of parameter space, c u ~ 1/2 (composite t/^), the situation is inverted, with the first KK 
modes of the members of the bidoublet of £~ becoming light, while the singlet is considerably 
heavier. These Z2-odd quarks couple to X_ and SM fermions (top or bottom), and we will 
refer to the lightest of them as the NLOP (Next-to-Lightest Odd Particle). Depending on the 
degree of degeneracy, coannihilation effects can be relevant to obtain the DM relic abundance. 
Note that for c u > 0, there are several states whose masses are split only by EWSB effects, 
and all of them can affect the final DM relic density. 

In order to get an order of magnitude estimate, we start by assuming that coannihilation 
processes can be neglected (a reasonable approximation whenever the Z2-odd quarks are 
heavier than X_ by a factor of > 15%). In this case, the only relevant process for the 
calculation of the DM relic abundance is the annihilation of X- pairs into SM quarks, via 
NLOP exchange. Up to EWSB effects, which can be shown to give negligible corrections, the 
relevant 5D interactions are 

(M)^-(l,l)^+h.c, (3.4) 



for c u ~ —1/2, and 



(2,2)^_(2,2)^+h.c., (3-5) 



for c u ~ 1/2. Here (1,1)^ and (2,2)^ contain, respectively, the right-handed (RH) top 
and the left-handed (LH) top and bottom zero modes. We denote by gn and g~L the 4D 
couplings between the LOP, NLOP and top/bottom quarks, arising from Eqs. (3.4) and 
(3.5), respectively. These are the crucial couplings entering in the self-annihilation cross 
section of X_. In the absence of UV brane kinetic terms, both gi and gn are bound from 
above as follows: 

/ 2 2 gg' , . 

\~9l,r\ < ^gx = r . = = ~ 0.28 , (3.6) 

where gx is the U{l)x gauge coupling that we have also written in terms of the SM SU (2) l x 
U(1)y gauge couplings. We refer the reader to appendix [C] for a detailed derivation of the 
couplings gi and c/r. They are typically of the same size, so that the main difference between 
the above two regions is the presence of one decay channel (into i/j) for c u < 0, and two (into 
ti and bi) in the c u > case. For this reason, it will be useful in the following to define the 
parameter rj, such that rj = 1 for c u < and r\ = 2 for c u > 0, and write the coupling simply 
as g, with the understanding that g = (Jr for c u < and g = gi for c u > 0. 

The computation of the DM relic abundance is standard (see e.g. [33] )■ In the freeze-out 
approximation, it can be written in terms of the coefficients of the non-relativistic expansion 
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of the annihilation cross section, <r(X_X_ — * SM), 

va = a + v 2 b + ..., (3.7) 

as 

Afp y^a + 3b/x F ( ' 

where Mp ~ 1.22 x 10 19 GeV is the Planck mass, xp = mx_/Tp, with Tp the freeze-out 
temperature (in our case, with < TeV weakly interacting particles, xp ~ 24 — 26), <?* is 
the effective number of relativistic degrees of freedom at freeze-out (g* = 86.25 for 100 GeV 
i$ i$ TeV) and a, 6 are measured in GeV -2 . Neglecting EWSB effects, the annihilation 
cross section depends on g, mx_ and the mass of the lightest modes in £~, which we call mwj. 
It turns out that the factor b gives a negligible effect. Denoting 

A = ^" m ^ , (3.9) 



we obtain 



and 



2r?g 4 1 _ 2 

3vr (1 + (1 + A) 2 )2 mx - ' 

\2 \ 2 



(3.10) 



2 1X15 / m x _ y /028V / 1 + (1 + A)2 \ 
77 UoOGeVv 1 V 9 J \1 + (1 + 0.15)2 J " 
This result is actually very accurate. Using the full expression for the annihilation cross 
section and the freeze-out temperature, we obtain, again assuming that the difference in 
mass between X_ and the Z2-odd quarks is > 15%, 

, f 0.16 (0.6), C „<0, 
[ 0.08 (0.28), c u > , 

where we assumed mj_ > 400 GeV, and the numbers are for g = 0.28 (maximal coupling) 
and, in parenthesis, for g = 0.2. We therefore see that the observed relic abundance [35 , 

n DM h 2 = 0.1143 ± 0.0034, (3.13) 

can be accounted for with just X_ annihilation for c u > 0, provided g is large enough. For 
c u < 0, the observed DM energy density seems more difficult to accommodate in the simplest 



scenario where only V_ annihilations are relevant 10 However, as explained above, EWSB 

10 We note that the presence of kinetic mixing between the two U(l) factors, or UV brane localized kinetic 
terms for them, can change the coupling relevant for annihilation of LOP's, leading to a realistic DM relic 
abundance. We assume that such terms are small, so that no relevant additional parameters are introduced. 
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and the EW precision measurements independently point to a region in parameter space 
with one or more vector-like quarks that are nearly degenerate with the LOP, and a proper 
determination of the relic abundance should take coannihilation effects into account [36]. 
This is welcome because, due to the colored nature of the new particles, coannihilations tend 
to increase the cross section and therefore decrease the relic abundance to values compatible 
with observation for smaller couplings and/or larger mx_- 

The relevant processes that enter the relic abundance computation are X_X_ — > ti(bb), 
X_ip_ — > gt (gb), X-tjj- — > gt(gb), ip_ip_ —* qq and if)_ip_ — > gg, where ip- stands for any of 
the Z2-odd fermions that are nearly degenerate with X_, g is the SM gluon and qq are SM 
quark- antiquark pairs. The relevant effective annihilation cross section reads 



1 



VcE — ~2~ 
9 C R 



+ 9(r? - 1) <? x x _^ b i 
+ 72 e~* A (l + Af/ 2 a x _^ gt + 72( V - 1) e^ A (l + Af/ 2 a x „ 
+ 72^-^(1 + A) 3 (j> v 

g cS = 3 + 12 7^-^(1 + A) 3 / 2 , 



(3.14) 



where x = mx_/T, and A was defined in Eq. (3.9). In Eq. (3.14), the factor rj introduced 



before takes into account that in the c u < region there is a single Z2-odd fermion with 
~ mx_, whereas for c u > there are four fermions with ~ mx_- When c u > 0, 



Eq. (3.14) is valid in the limit in which we neglect the mass splitting between the fermions 



of the bidoublet (which arises only from EWSB effects) and it is understood that the cross 
sections with final state top (bottom) quarks involve the Q = 2/3 {Q = —1/3) heavy vector- 
like fermions. 

We performed exhaustive scans, computing the pattern of EWSB, the fit to EW pre- 
cision data, the DM relic abundance and the spectrum of light quarks. In order to more 
easily automatize the scan over large regions of parameter space while taking into account 
coannihilations, we found it useful to implement the relevant features of our model in mi- 
crOmegas [37 . This also allows us to easily take into account the mass splittings arising from 
EWSB among the lightest bidoublet states that are relevant in the c u ~ 1/2 region, as well 
as non-zero final state masses (for the typical scales of mj_ ~ 400 GeV, the top mass can 
be important). We have checked independently that we reproduce the micrOmegas results 



for the relic density to within 10% using Eq. (3.14), in the limit that EWSB mass splitting 



effects are neglected, and in the freeze-out approximation. 

We present the results separately for the two phenomenologically interesting regions dis- 
cussed previously, namely c u ~ —1/2 and c u ~ 1/2. We scanned over the following region in 
parameter space: c qi & [0,0.4], c q2 G [0.4,0.5], \c u \ € [0.45,0.5], Cd G [-0.52,-0.45], \m u \ < 1, 
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Figure 1: Projection onto the c u -c qi plane for the two phenomenologically viable regions. 
The darker areas marked Sh = correspond to no EWSB. We plot Sh averaged over the 
rest of the parameters, which increases as the gray bands become lighter (see text). The 
green dots satisfy 145 GeV < mt(/i ~ Mr) < 155 GeV. The red triangles mark the points 
consistent with the WMAP constraint, Eq. ( |3.13 ), at the 2a level. The blue stars correspond 
to a sample of points that are consistent with EW precision data at the 99% CL, and the 
Higgs LEP bound. 



\M U \ < 1, |m<i| < 0.5, \Md\ < 0.5, with kL ~ 34 in order to explain the Planck-weak scale 
hierarchy. We show the results of the sccin as 3. projection onto the c u -Cq^ plane in Fig. [TJ 
where the red triangles correspond to the data consistent with the three- year WMAP results, 
Eq. pl3| ), at the 2a level. 

We also highlight in the plot the data that have 145 GeV < m t (fi ~ /i m ) < 155 GeV (green 
dots). As we will discuss further in section 3.4, the top mass plays a rather important role in 
leading to a vacuum with appropriate characteristics in these scenarios. In the background 
of the figure we exhibit in gray tones information about the degree of EWSB as measured by 
Sh [see Eq. (3.1)]. In order to make the projection onto the c u -c qi plane we show the average 
of Sh, computed over the rest of the parameters in the scan. In particular, in the darker area 
marked as Sh = 0, all points in the scan lead to an EW symmetry preserving vacuum. For 
c u < (c u > 0), the degree of EWSB, in the above sense, increases as one moves towards the 
bottom left (right) corner of the figure (lighter grays correspond to larger average Sh)- We 
therefore see that in most of the region with non-trivial EWSB the correct top mass can be 
reproduced. Even more interestingly, the area selected by the observed DM abundance falls 
in precisely the same region. This is related to the fact that coannhilations play a relevant 
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Figure 2: Masses of the DM candidate, X-, and the Higgs for points (red dots) that reproduce 
the top mass and the WMAP constraint, Eq. (3.13), for the c u < region (left panel) and 
c u > region (right panel) . The blue stars correspond to the subset of points that also obey 
the EW precision constraints at the 99% CL, and the LEP bound on the Higgs mass, which 
is also indicated by the dashed horizontal line. The light bands indicate the approximate 
bounds from the Tevatron on the colored vector-like quarks, that have mass close to mx_, 
taking into account the different multiplicities for positive and negative c u . 



role in lowering the DM abundance to the observed level compared to the estimates Eq. (3.12 ) 
and, as explained before, this happens naturally in the above region (we did not compute the 
DM relic abundance for points with Sh = 0). 

Finally, we also show in the figure information about the fit to the EW precision data. As 
explained in section 3.1, we include in our fit the universal corrections described by the S and 
T parameters and the correction to the Zb^bi vertex (both tree- and one-loop contributions), 
performing a A% 2 cut at 99% CL with 3 d.o.f. Since this analysis is computationally intensive, 
we performed the EW test only for points that satisfy 149.5 GeV < mt(/i ~ fijR.) < 150.5 GeV 
(without imposing the WMAP constraint), and for 145 GeV < mj(/i ~ (J>ir) < 155 GeV when 
the WMAP constraint, Eq. (3.13), is imposed. The selection of the very narrow range in the 
first case is done only in order to reduce the number of points to be analyzed, but it should 
nevertheless give a clear picture of the situation. We also require that the Higgs mass be 
above the LEP bound, nih > 114.4 GeV. Requiring that the DM relic abundance be in the 
range (3.13), we find that in the c u > (c u < 0) region about 25% (15%) of the points pass 
the EW precision test at the 99% CL. Thus, we conclude that in the present model, the 
correct EWSB pattern, consistency with the EW precision measurements, and the correct 
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DM abundance all occur in a common region of parameter space. This leads to a rather 
compelling picture with a light neutral spin-1 DM candidate, and one or more fermionic 
resonances that are nearly degenerate with the DM particle. We will briefly discuss their 
collider phenomenology in section [5] 

Two important physical observables that are relatively well predicted in our scenario are 
the Higgs mass and the mass of the DM candidate. We show in Fig. [2] the distribution of 
points with a realistic EWSB pattern (including the top mass) that reproduce the observed 
DM relic density, in the mx_ — rrih plane. We also highlight the subset of points that are 
consistent with the EW precision constraints at the 99% CL. We see that the LEP bound 
on the Higgs mass is easily evaded in our scenario. Also, the Higgs boson is expected to 
be lighter than about 170 GeV. The mass of the DM particle is expected to be somewhat 
above 300 GeV (400 GeV) for c u > (c u < 0). As was mentioned above, there is a number of 
vector-like quarks nearly degenerate with the LOP. We also indicate in Fig.[2]the approximate 
direct Tevatron bound on such particles. Although the direct bound already excludes a 
significant region in parameter space, the EW precision analysis still puts stronger constraints 
on the masses of these particles. It is also useful to remember that, for kL ~ 34, we have 
mj_ 0.24 fijji, which shows that /j,jr > 1.3 TeV ((J-ir > 1.7 TeV) for c u > (c u < 0). 
Recall also that, in the absence of brane kinetic terms, other spin-1 resonances have a mass 
m KK « 2.5 mr > 3.1 TeV (m KK > 4.2 TeV) for c u > (c u < 0). 

3.3 Non-perturbative Corrections 

There are various physical processes that might possibly invalidate the standard computation 
of the relic density based on perturbative averaged cross-sections, namely the formation of 
bound states and higher order corrections relevant for non-relativistic particles. Such effects 
do not directly affect the physics of our DM candidate X_ , but they can affect the physics of 
the NLOP, which is a charged and colored fermion. Since coannihilations have to be taken 
into account in our set-up, it is important to estimate the above effects. We will not perform 
a detailed quantitative study of bound state formation at finite temperature since, as we 
will argue below based on qualitative estimates, such effects can most likely be neglected. 
On the other hand, certain higher order corrections, which take into account the long-range 
Coulomb-like forces of QCD in the deconfinement phase, are potentially relevant. We find, 
however, that their eventual contribution to the perturbative cross-sections is very small. 
Since strong interactions play the major role, we will neglect in this subsection the effects 
due to electroweak interactions. 

Let us first consider the issues associated with the formation of bound states. The poten- 
tial problem is very simple. If meta-stable, Z2-even, bound states of two NLOP particles can 
form, their decay through the self-annihilation of its constituents would lead to an effective 
depletion of Z2-odd particles, consequently reducing the DM relic density. Bound state ef- 
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fects of DM colored particles (gluinos) have been previously considered in [38 , in connection 
with relic density computations, and were shown to have the potential to reduce the final 
DM density by orders of magnitude compared to the perturbative estimate. The crucial dif- 
ference here is that the NLOP's are unstable, and as we argue below they decay into LOP's 
well before they have time to form bound states. The NLOP can decay via inter-generational 
mixing into the LOP and a light quark (decays into tops are forbidden by phase space in the 



quasi-degenerate region of interest to us). To lowest order in A, as defined in Eq. (3.9), the 
NLOP lifetime is 

Snm- 1 2 x 10~ 23 

T ^ 3A2^^ A 2 S6C ' ( ^ 

where A 2 = J2 q =cu(.9qL + 9qn) 1S the (model-dependent) coupling of the NLOP to the up 
and charm quarks, and the number quoted is obtained by taking a rather degenerate case, 
mx_ = 350 GeV and = 360 GeV. 

We estimate next the time scale when the (NLOP-NLOP) bound states would become 
meta-stable. A reasonable criterion is to assume that these bound states are meta-stable 
against thermal fluctuations when the temperature drops below their binding energy E^nd- 
Since the bound state system is non-relativistic, and in fact the scales are such that the QCD 
interactions are in the perturbative regime, we can estimate Ef^nd in analogy to positron- 
ium. [^] The color Coulomb force between two quarks in a color singlet state are obtained 
from the usual Coulomb potential by the replacement a — > (4/3) a s (p). Here a s (p) is eval- 
uated at the typical momentum scale of the virtual gluons responsible for the Coulomb-like 
interactions, which is p ~ m^a s (p)/2. For ~ 400 GeV -the typical scale in our scenario- 
one gets p ~ 30 GeV and ce s (p) — 0.14. The binding energy is then given by 

Ebind - ^ (^) a s (p) 2 m^ ~ 10 _2 m^ . (3.16) 

Thus, the (NLOP-NLOP) bound states become meta-stable after freeze out, which occurs 
when T ~ m^/25. In this radiation dominated era, these two events happen at times t ~ 10~ 7 
sec. and t ~ 4 x 10~ 9 sec. in the evolution of the universe, respectively. It is therefore 



clear from Eq. (3.15) that for any reasonable value of the model-dependent coupling A, the 
NLOP's decay into LOP's immediately after freeze-out, and well before bound states can 
become meta-stable. 

The second effect we consider is due to the long-range Coulomb interactions that, for 
sufficiently non-relativistic particles, can distort their wave functions from the plane wave 
shape and alter the standard quantum field theory perturbative computation of scattering 
processes. Such an effect is well-known and was first analyzed by Sommerfeld in the QED 



"Strictly speaking, at finite temperature there is a screening effect on charged particles, leading to an 
effective Debye mass for the gluon, m g ~ V2g a (T)T. It is easy to check, however, that for T/m^ ~ fO [see 



Eq. ( 3.16 )] , the thermal screening can be neglected. 
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context [39]. It has recently been considered, in connection with relic density computations, 
in jlQ] for electroweak interactions and in [38] for strong interactions. For QED at zero 
temperature, the "Sommerfeld" effect can be encoded in an effective parameter S which 
reads, for absolutely stable particles P| 



S 



±x 



1 



3 ±X 



(3.17) 



with x = ira/v, v the velocity of the colliding particles in the center of mass frame, and ± 
refer to repulsive or attractive Coulomb forces, respectively. Given a cross-section a com- 
puted in the standard perturbative fashion between charged non-relativistic initial states, the 
replacement a — * Sa is an effective way to take into account the Coulomb forces, assuming 
relativistic final states. In a relativistic approach, the Sommerfeld factor S is obtained by 



resumming an infinite class of Feynman diagrams ("ladder" diagrams). From Eq. (3.17), it is 
clear that S is non-negligible only for sufficiently non-relativistic particles. Cold DM candi- 
dates are by definition non-relativistic at freeze-out and the Sommerfeld effect can play a role, 
as emphasized in |38| l4"0] . For strong interactions in a perturbative regime, an analysis along 
the lines of QED can be made. By neglecting screening effects due to temperature (which is a 



good approximation at Tf, see also footnote 11), Eq. (3.17) still holds, with the replacement 
a — ► C r a s (p) in x, where C r is a color factor that depends on the SU(3) representation of 
the particles involved, and p ~ mv is the typical momentum of the virtual gluons responsible 
for the Sommerfeld effect in a relativistic treatment |38J . For a tp-ijj pair in the initial state, 
one can have either a singlet or an octet SU(3) configuration, with color factors C\ = 4/3 
and Cg = —1/6, respectively At freeze-out, v ~ ^2Tf /m^ ~ 0.25 and we get, for ~ 400 
GeV, Si — 2.3 and Sg — 0.9. Thus, the Sommerfeld effect, especially in the singlet channel, 
is non-negligible. We have quantified the impact of the Sommerfeld effect on the computation 
of £1 by replacing the averaged tree-level cross-sections as follows: 



(7 



i/jjp—t gluons 



a 



ipi/j^qq 



uons ^ ^8 l -f(ip'tp') 8 —> gluons 
- + 8 S 8 V(W) a ^qq 



As it turns out, the enhancement of the cross-sections in the attractive singlet channel is 
largely compensated by the repulsive octet channel, so that the effect on the final O is always 
very small (at most a few percent). 



12 For unstable particles, one has to check whether the Sommerfeld effect has time to take place, by comparing 
the typical time scale of Coulomb interactions tcouiomb ~ l/(Mv 2 ) with the decay time r of the particle |41] , 
If t <C tcouiomb, the Sommerfeld effect has no time to take place. In our case, tnlop 3> tcouiomb and thus 
we can effectively treat the initial NLOP particles as absolutely stable. 
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3.4 Fine-tuning 

We established above that a rather well-defined region with \c u \ ~ 1/2 is simultaneously 
selected by the minimization of the Higgs potential, the EW constraints and the observed 
DM relic density. We should also recall that the pseudo-Goldstone nature of the Higgs allows 
for its mass to be parametrically lower than the scale of the KK resonances /i IR , potentially 
alleviating the little hierarchy problem present in other RS constructions (of course, the large 
Planck- weak scale hierarchy is explained by the RS mechanism). It is therefore natural to 
ask how fine-tuned these scenarios really are. We quantify the fine-tuning at a given point in 
parameter space by considering the sensitivity to the microscopic parameters of the theory 
as measured by the logarithmic derivative [H] 



sensitivity = max 



dlog(h) 
d log Xi 



(3.18) 



where Aj = c qi , c q2 ,c u ,Cd, m u , M u , ma, are the fundamental parameters of the model. We 
find that the apparent fine-tuning is dominated by the sensitivity to c u , as expected from our 
previous discussions. We show in the left panel of Fig. [3] the sensitivity parameter for the 



random scans described in section 3.2 as a function of Sh- These show that throughout the 
region of parameter space that leads to phenomenologically viable EWSB breaking minima, 
the apparent fine-tuning is worse than a percent. Without tractable analytic expressions 
it is hard to identify precisely the nature of such sensitivity, and it is conceivable that the 
exponential nature of the warp factor gives rise to a generic sensitivity. In a 4D dual language, 
it would probably be related to the usual generic sensitivity which affects theories with a 
dynamical generation of scales, such as QCD or technicolor theories |43j . It is then possible 



that the fine-tuning estimate as given by Eq. (3.18) is too conservative and a more refined 
analysis is necessary. It would be interesting to investigate in more detail this issue. 

However, we notice here that the sensitivity to the microscopic parameters seems to 
be almost exclusively related to the requirement of obtaining a large enough top Yukawa 
coupling. In order to illustrate the point, we compute again the logarithmic derivatives in 
the (c u , M u ) plane, along directions with constant top mass. The result is shown in the right 
panel of Fig. [3j Although we have not explored the logarithmic derivatives at constant top 
mass in the whole parameter space and thus cannot exclude the existence of other directions 
where the sensitivity is higher, Fig. [3] shows that part of the fine-tuning in the model is 
associated with reproducing the observed top mass. Similar results also hold in the original 
MCHM5 model without DM. We therefore conclude that the present model necessitates some 
degree of fine-tuning to reproduce the observed top mass, but once this measurement has been 
done, the physical properties of the resulting vacuum are fairly robust against variations in the 
microscopic parameters of the theory (observables other than (h) also exhibit this property). 
Notice that this is a single measurement associated with the zero-mode sector, and therefore 
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Sh s h 



Figure 3: Left panel: logarithmic sensitivity as a function of s^. Right panel: logarithmic 
sensitivity with the requirement that the top mass be kept fixed (see text). The green dots 
correspond to a random scan over parameter space with c u > 0. The subset that satisfies 
140 GeV < mt(/j ~ fim) < 160 GeV is indicated by red triangles. The results for the c u < 
region are similar. 

there is a sense in which there is little sensitivity to the details of the new physics beyond the 
standard model. This also underscores the role that the heavy top plays in driving EWSB. 

The above situation should be contrasted with supersymmetric (SUSY) scenarios that 
also present, generically, a sensitivity of order a percent to various fundamental parameters. 
In the SUSY case, reproducing a large top mass is no more difficult than in the SM (it is 
enough to choose the top Yukawa coupling, one of the "fundamental" parameters of these 
theories, to be of order one). However, there is a well-known intrinsic fine-tuning associated 
with a cancellation between the \i term and the soft SUSY breaking parameters in the Higgs 
sector. The latter depend quadratically on the stop soft SUSY breaking masses, which in turn 
need to be taken somewhat heavy (at least in the simplest SUSY extensions of the standard 
model) in order to get a large enough Higgs quartic coupling that allows satisfying the LEP 
bound on the Higgs mass. If, in addition, one takes into account the RG running from a high 
scale, and quantifies the fine-tuning by the sensitivity to the high- energy parameters of the 
theory, the situation worsens. In warped scenarios, the RS mechanism eliminates any possible 
fine-tuning due to running from a high scale. The dynamical generation of the weak scale in 
gauge-Higgs unification scenarios, allows for a further natural separation between the weak 
scale and the KK scale fi m . In fact, EWSB is generic, except that in most of parameter space 
it is characterized by = 1, which is not phenomenologically acceptable. The sensitivity 
shown in Fig. [3] is associated with the requirement < a/, < 1. Our observation is that all 
of this sensitivity is actually associated with getting the correct top mass (which in these 
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scenarios is related in a relatively complicated manner to the fundamental parameters of the 
model), rather than with the LEP bound on the Higgs mass, which as Fig. [2] shows is easily 
satisfied in these scenarios. 



4 Dark Matter in Other Models with Warped Extra Dimensions 

In the previous section, we showed how a very simple extension of a minimal composite Higgs 
model in warped extra dimensions successfully accounts for the observed DM relic abundance, 
while leading to a non-trivial connection with the physics of EWSB and EW constraints. We 
now briefly show that our prescription can be easily implemented in essentially any model 
with warped extra dimensions. 



4.1 RS Model with Bulk Gauge and Fermion Fields 

This is the simplest model in which our mechanism can be applied. We will see that it is 
possible to get the correct DM relic density, although some fine-tuning is likely required. The 
bulk gauge symmetry is taken to be SU(3) C x SU(2)l X U(l)x 1 x U(l)x 2 with the Higgs field 
localized on or near the IR brane. The SM U (1) hypercharge factor is identified with the even 
linear combination, U(l)y = U(l)x + - The DM candidate X- is the first KK mode of the 5D 



odd combination, Eq. ( |2.4| , and has a mass given by Eq. (2.5 ), ~ y/2/kL fj, m ~ 0.24 /j, m , 
where we have used kL ~ log (Mp /TeV) ~ 34. Assuming no brane kinetic terms, the IR scale 
in this model must obey // IR > 3 TeV |2j in order to be consistent with the EW precision 
measurements. 13 The annihilation rate of X_ is again fixed by the 4D coupling g between 



X-, the lightest Z2-odd fermion, ip~ \ an d the SM zero-mode field ip+, an d their masses. 
A simple computation shows that g < Yg', with g' and Y the SM hypercharge coupling and 
hypercharge quantum number of ip+^- The Z2-odd fermions arise as explained in section [2] 
by doubling a subset of the (SM) fermion fields. These must be chosen to have a localization 
parameter |c| < 1/2, so that the resulting Z2-odd fields are heavier than JT_. In scenarios that 
explain the fermion mass hierarchies by fermion localization, this leaves as natural candidates 
the quarks in the third generation. The self-annihilation rate cr(X_X_ — ► is too 

small to give the correct relic density f2/i 2 , which is correspondingly too high. However, the 
DM relic density can be lowered to the observed level provided coannihilation processes with 
the colored quarks tp^ are important. This requires ip_ and X- to be fairly degenerate, 
hence |c| to be very close (from below) to 1/2. It turns out to be hard to take c qL ~ 1/2 
or ct R — —1/2, because in order to get the correct 4D top Yukawa coupling one is forced 
to increase so much the 5D top Yukawa coupling that it would enter the strong coupling 



13 Allowing for the Higgs to be a bulk field with an exponential localization towards the IR brane [45] 
reduces this lower bound to fj, m > 2 TeV. Similarly, sizable IR brane kinetic terms for the gauge bosons can 
also decrease the corresponding bound on the masses of the KK modes significantly |46] . 
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regime [3 . The alternative is to double the right-handed bottom quark only, and take q> H — 
—1/2. In this case, other scattering processes become important, namely ct(X_&r_ — > b^g), 
a{bR-bR- — > gg) and <t(6r_&r_ — > qq). Given enough degeneracy between X_ and the 
relic density of X_ can match the observation. For instance, for mx_ = 750 (500) GeV, 
which corresponds to /Utr. — 3 (2) TeV (see footnote 13), one requires ci R ^ —0.496 (—0.495). 
In terms of mass splittings, one gets A = — mx_)/mx_ — 6 (9)%. This is not a 

fully satisfactory scenario due to the moderate tuning needed (and the lack of independent 
motivation for the choice of parameters), but it shows how this idea can be successfully 
implemented even in the simplest 5D warped model. 

By relaxing the explanation of hierarchical Yukawa couplings due to the fermion localiza- 
tion one can envisage a different scenario, where all (light) SM fermions share the same profile 
in the extra dimension. Such a situation occurs in Higgsless models, which is the subject of 
the next subsection. 



4.2 Higgsless Models 

Higgsless models induce EWSB by means of boundary conditions in the extra dimension [U] . 
In this way, the SM gauge bosons acquire their longitudinal components through the Higgs 
mechanism but no extra scalar degree of freedom (the Higgs) is present in the low energy 
spectrum. Unitarity violations in longitudinal gauge boson scattering are delayed by the 
exchange of the gauge boson KK excitations up to the cut-off of the theory. Recently, using 
ideas borrowed from [3j 3] , the first 5-dimensional Higgsless model roughly compatible with 
EW precision data at tree level has been presented jH]. [^] A crucial ingredient is the 
derealization of the (LH) light fermions [50 , that have a common localization parameter 
close to the conformal point cl < 1/2. This means that we can double the whole (LH) 
fermionic spectrum, following the prescription described in section [2j with the corresponding 
increase in the number of open channels for the DM particle to annihilate. Furthermore, 
the low IR scale favors a stronger annihilation cross-section and naturally produces the right 
order of magnitude for the relic abundance. Our starting point is the model of Ref. |48] . 
in which cl = 0.46. In this case, the bulk gauge group is SU(2)l x SU(2)r x U{l)x^ x 
U(1)x 2 i with U(1)b-l = U(l)x + the even combination of X\ and X2. The smaller warp 
factor and lower IR scale in this model lead to a relatively light DM candidate, with mass 
mx_ ~ 0.4 /x IR ~ 114 GeV. Doubling the whole LH light spectrum leads to a too large 
annihilation cross section, mostly due to the larger charge of the leptonic fields under the 
U(1)b-l g rou P) whereas doubling only the light LH quarks results in a cross-section that is 
a bit too small. As an example, doubling only one LH lepton doublet, we obtain £lh 2 ~ 0.07. 

14 Deconstructed versions of Higgsless models are more flexible than 5-dimensional ones and models com- 
patible with data, at the tree level, have been constructed [49j. 
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This is impressively close to the observed value, considering that we have taken identical 
values for all the parameters, as in Ref. [48 . 



5 Phenomenology of the Dark Matter Sector 
5.1 Collider Phenomenology 

The collider signatures of our construction share common qualitative features with other 
models of new physics with stable particles. New states are produced in pairs and (cascade) 
decay to SM particles and the LOP, which is perceived as missing energy in the detector. 
The specific details of the spectrum of NLOP induce, however, significant differences at the 
quantitative level. 

In order to be more specific, we will discuss the phenomenology of the GHU model 
described in section [3j We study, as an interesting and representative example, the case 
with c u > 0. The light Z2— odd spectrum consists of the LOP with a mass mx_ ~ 300 — 400 
GeV and the first KK modes of the bidoublet component of £~, that contain two quarks 
of charge 2/3, one of charge 5/3 and one of charge —1/3. We will denote these quarks as 
t', q u , x u an d q d , respectively, and, generically, as ip. They are degenerate, with a mass 
m^p — mx_ ~ 0.15 mx_, except for small EWSB effects that make one of the charge 2/3 
quarks (t') slightly lighter. 

Vector-like quarks are pair-produced at the LHC via the QCD interactions with a model- 
independent cross-section that depends only on the mass of the quark [51] , 

°pp->QQ ~ 15 (1.5) pb, for m Q ~ 400 (600) GeV. (5.1) 

The distinctive feature is the high degree of degeneracy between the NLOP and the LOP, 
which forbids the natural cascade decay through a top quark, leaving decays into light jets 
and missing energy (the LOP) as the main signature. The charge 5/3 quark cannot directly 
decay into the LOP and SM particles, and therefore undertakes a cascade decay through an 
off-shell W* to t', which then decays to jets and missing energy. 

The decay width for the two-body decays (kinematically allowed for q d and, through 
inter-generational mixing for t' and q u ) is given, to leading order in the mass difference, by 



(the inverse of) Eq. (3.15) 



T U - j X J) « — A 2 m^ «(4x 10~ 6 GeV) ( —) ' ( —) ' (5.2) 
^ J ' 8vr * v ; V - 002 / V - 15 / 460 

where j = u,c,b, the relevant coupling is denoted with A, ip stands here for any of t',q u 
and q d , and is measured in GeV. For the numerical result we have replaced the different 
parameters with typical values for the charge 2/3 quarks (we have assumed inter-generational 
mixing in the up sector to be of the order of the corresponding Cabbibo-Kobayashi-Maskawa 
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matrix elements A ~ V u b iC b9/V^ ^ 0.2 x (10 2 — 10 3 )). In the case of q d there is no 
inter-generational mixing suppression (A = g) and the width can easily be 

T{q d bX_) » (lO^ 1 - 10~ 2 ) GeV. (5.3) 

The corresponding three-body decay for \ u proceeds through an off-shell W* , which can 
decay leptonically. In the contact interaction approximation, the partial decay width into 
muons is 

nr - G -Mf m « ( . * *>- cv, <m> 

where f(x) = 1 — 8x 2 + 8x 6 — x 8 — 12x 4 logx 2 , with x = mf /my*. 

Thus, we see that the NLOP are produced with a large cross section and decay promptly, 
leading to the following typical signatures, 

pp - t'V, q u q», q d q d - jj fa, (5.5) 
- X a X u - - Zi/jjij #r- (5.6) 

The former signature, that can benefit from a larger cross-section due to the quark multi- 
plicity, is challenging due to the lack of leptons to trigger on, and the fact that the amount 
of transverse missing energy is limited by the small mass difference between the NLOP and 
the LOP. The latter, more promising due to the presence of a leptonically decaying W, has 
however the problem of the extra source of missing energy (the neutrino from the W*) and 
the fact that the W* is off its mass shell and therefore its mass cannot be reconstructed. 
A detailed analysis, which is beyond the scope of this work, is needed to asses the trigger 
efficiency in accepting the signal and the best strategy to search for this new sector in models 
with warped extra dimensions, but the large production cross-sections seem to indicate that 
discovery should be possible at the LHC. 

Thus, although the generic collider implications of our construction share some features 
with other models of new physics with stable particles, the particular details of the spectrum 
of NLOP gives a very characteristic signature, with direct decays into jets plus missing energy 
and, in some cases, also short cascade decays through off-shell Ws. These are challenging 
signals at the LHC, but the large cross sections due to the low mass of the new particles 
and the high multiplicity should help in the discovery of these channels. One should also 
remember that the Z2-even sector of these models is also predicted to have light resonances. 
In the particular case we are discussing, we have a light quark doublet with hypercharge 
7/6 that will be easily observed in the early phase of the LHC [52] . The more challenging, 
although possible [53 j , discovery of the gauge boson KK excitations and the lightest particles 
of the Z2-odd sector as we have described above, should then draw a clear picture of the 
structure of EWSB. 
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5.2 Direct Dark Matter Searches 



Direct detection of DM particles is mostly based on detectors in which the relevant process 
is the scattering amplitude DM-atoms^ DM-atoms. Since, as we have seen, X_ typically 
couples only to third generation quarks (the top quark in the GHU model, bn in the RS 
model), its direct detection rate is expected to be too small for current and future planned 
experiments. We can be more quantitative and estimate the cross-section X-N — * X-N, 
where 7V is a nucleon. Not surprisingly, the situation is analogous to that of UED where the 
DM particle is identified with the first KK mode of the 5D hypercharge gauge field, with the 
important difference that in our case there is no Higgs exchange and only heavy quarks Q are 
relevant. Following [53] and using their notation, the spin-independent cross-section reads 

tit? 

^scalar = "j 7 . vj /jV > (5-7) 

47r(mx_ + mN) 
where Jm = /3q(N\QQ\N) with Q top or bottom quarks, and 

9 9 
TJX ~\~ 7Tb 

0Q - m Q (Q 2 XjL g 2 L + Qx,r9r) , 2 t 2 , (5-8) 

\m x _ m^) 

with m^p the mass of the NLOP. A careful estimate of the nuclear matrix element (N\QQ\N) 
would require a detailed one-loop analysis, along the lines of |55] . For an order of magnitude 
estimate, it is however enough to use the old result [56] 

(mQQ\N) = ^(l- £ /£), (5.9) 

q=u,d,s 

where = {N\qq\N)m q /mj<i . It is now straightforward, to compute o scalar • 

For the GHU 

model we considered, by taking, say, mx- — 350 GeV, ~ 380 GeV, = 0, gn ~ 0.25, 
Qx,R = 2/3, one has 

O scalar -2 X 1(T 10 pb , (5.10) 

which is a value too low for current experiments and would require new experiment proposals, 
such as super-CDMS [57] . Due to the higher value of mj_ , a sca i ar is even smaller in the RS 
model with doubled bn considered in subsection |4.1[ In both models, due to the denominator 



term in Eq. (5.8), cr scalar can become sizable only for extreme degenerate cases in which 
m x _ ^ m^j. 

In Higgsless model, the situation is different, since one might also have direct couplings 
of X- with light quarks. In addition, X- is lighter, m x _ ~ 100 GeV, and hence DM 
direct detection S66ins more promising. However, in the specific model of |48j, cr scalar is still 
suppressed due to the low values of the charges, Qx,L = 1/6 (Qx,R = 0), resulting in a 



cross-section of the same order of magnitude as (5.10) 
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6 A Comment on Anomalies 



We have assumed so far that the Z2 exchange symmetry is an exact symmetry of the theory, 
namely that no quantum corrections can possibly violate it. We show here that the exchange 
symmetry is exact by noticing that the CS terms that are required by gauge invariance are 
always Z2-even and hence invariant. This is actually expected, since the Z2 symmetry is a 
global symmetry which has nothing to do with parity, the discrete symmetry typically broken 
by anomalies. For simplicity, we will focus our attention on the simplest RS scenario with 
fermion and gauge bulk fields analyzed in subsection 4.1 but the main results remain valid 
also for the more refined GHU model of section |3j 

The bulk gauge group is G = SU{3) C x SU(2) L x U(l) Xl x U(l) Xa> broken to SU(3) C x 
SU{2)l x U(l)y at the IR brane and fully unbroken at the UV brane. The 5D fermion 

spectrum consists of massive Dirac fermions, one for each SM fermion, with (++) or ( ) 

b.c, depending on the chirality of the SM fermion. They have charge {\Y, ^Y) under U(\)x\ x 
U(l)x 2 i with Y the corresponding SM hypercharge for the given fermion field, while the 5D 
bottom fermion fields ip^ and ipb 2 have charges (—1/3,0) and (0,-1/3). Although the 4D 
massless fermion spectrum of the model is anomaly free by construction, being identical to 
the SM spectrum, the 5D theory needs CS terms to restore 5D gauge invariance fully, because 
localized anomalies (globally vanishing in 4D once integrated over the internal space) appear 



at the UV and IR branes (see |58| for a review) 15 Recall that the warping does not alter the 
localization pattern of anomalies, which is then as in flat space [59]. It is extremely useful 
to also recall that the form of a localized anomaly at a boundary is fully determined by the 
"effective" chiral spectrum which is found there by neglecting the b.c. at the other end-point, 
which is equivalent to sending the other boundary to infinity (see e.g. [58] for a derivation of 
this result). Let us consider the UV brane, where both U(l) factors are unbroken, and focus 
on the possible anomalies involving the gauge field All SM fermions are neutral under 
U{l)x-i with the exception of tpb 1 and ipb 2 , which have opposite charges —1/3 and +1/3. 
Correspondingly, all U(l) x _, U(l)x_U(l) Y , SU(3) c U(l)x^ and mixed f/(l)x_ -gravitational 
anomalies trivially vanish at the UV brane. At the IR brane, no anomalies involving U(l)x- 
can appear, since the gauge field X- vanishes there. It is not difficult to compute all other 
localized anomalies. One finds that the "doubling" procedure induces pure SU(3) C , U(1) Y 
and mixed SU(3) c U(1)y , J7(l)y-gravitational anomalies at the UV and IR brane, as well as 
a U(1)yU(1) x anomaly at the UV brane. 11 these anomalies are cancelled by suitable 5D CS 



terms, which involve either (the former) or 2 (the latter) gauge fields A_. 6 Consequently 



all CS terms are Z2-even, as anticipated. A similar result is found for the GHU model of 



In this simple example, localized anomalies appear only due to the doubling we have performed. The 
original model is fully anomaly free since the SM anomaly cancellation applies also to the localized terms. 

16 Possible modifications to the wave functions of gauge fields induced by such CS terms can be safely 
neglected, being one-loop suppressed. 
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section |3j 



7 Discussion and Conclusions 

In this paper we have proposed a generic construction that allows to endow given models 
with a DM candidate. This is achieved through an extension in which the model acquires a 
Z2 exchange symmetry. The lightest Z2-odd particle is then absolutely stable. 

We have considered several models with warped extra dimensions. Although these scenar- 
ios are well-motivated extensions of the SM, explaining both the Planck-weak scale hierarchy 
as well as the flavor structure of the SM, they do not contain, generically, stable particles 
that can account for the observed DM component of the universe. We have shown that our 
mechanism can easily solve this deficiency. 

We have payed special attention to a class of warped scenarios that is particularly ap- 
pealing: GHU/composite Higgs scenarios. In this case, the Z2 structure responsible for the 
stability of the DM candidate is tightly connected to the physics that leads to the dynamical 
breaking of the EW symmetry. As in supersymmetry, the dark matter mass and couplings 
are intimately connected to the EW scale. In fact, not only does our construction not in- 
troduce new parameters, but there is a further sense in which it can be considered minimal. 
As was emphasized in the main text, the physics of EWSB in such models crucially depends 
on certain fields that have two properties: they are fermionic fields without zero-modes and, 
through their strong connection to the top sector, they give a significant -in fact, crucial- 
contribution to the Higgs potential. This is precisely the sector that gives rise to the Z2-odd 
particles that can lead to a realistic dark matter candidate. As a result, the Z2-odd sector, 
through its contribution to the Higgs potential, plays a key role in the realization of an EWSB 
vacuum with the desired physical properties. As explained in detail in the main text, the 
effects due to the Z2-odd sector go in the direction of alleviating the constraints imposed 
by low-energy measurements. The picture that emerges is rather compelling: the observed 
DM relic abundance and the very precise measurements of the EW observables at LEP and 
SLC point to a common region in parameter space. This region is characterized by several 
fermionic resonances nearly degenerate with the DM candidate (a spin-1 particle), which in 
turn is predicted to have a mass in the 300 — 500 GeV range. 

We have also pointed out that although the GHU scenarios (with or without DM) present 
a sensitivity of order one percent or worse with respect to microscopic parameters (larger 
than what naive considerations would indicate) , this sensitivity disappears once the top mass 
measurement is imposed. Thus, the moderate sensitivity of order a percent, that is present 
in virtually every extension of the SM that has any relation to the physics of EWSB acquires 
a new twist: it indicates an intrinsic difficulty in accommodating a heavy top. By contrast, 
once the top mass is fixed to a value of order the EW scale, the Higgs mass easily satisfies 
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the LEP bound, while typically lying below 170 GeV or so. 

Regarding the phenomenological implications of our construction, the low scale predicted 
for the new Z2— odd sector and, in particular, new vector-like quarks, guarantees large pro- 
duction cross-sections. Due to the particular features of the NLOP spectrum, which is very 
degenerate with the LOP, the most common signature will be jets plus missing energy accom- 
panied in some cases by the semileptonic decays of a pair of off-shell W* . These are challeng- 
ing signatures that will require dedicated analysis. The large production cross-sections and 
the information coming from the easier channels in the Z2— even sector should however be 
sufficient to guarantee a discovery of the DM sector at colliders. We have also seen that DM 
direct detection is not very promising in present or near future experiments, due to the fact 
that couplings to light valence quarks are suppressed by intergenerational mixing, leading to 
very small cross-sections. 
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A The MCHM 5 with Dark Matter 

In this section we summarize the model of GHU that was analyzed in the main text. The 
model has an SO(5) x U(l)xi x U(l)x 2 bulk gauge symmetry, with a discrete symmetry 
under which the two U(l) factors are exchanged. The Z2-even sector is identical to the 
MCHM5 of Ref. [25]. The SO(5)/ SO(4) directions are broken on both branes by choosing 
(— , — ) boundary conditions for the // components 

4 (-,-), (A.l) 

where a = 1,...,4 runs over the SO (5) /SO (4) indices. The corresponding components 
along the extra dimension, (which are four dimensional scalars) have zero modes that 
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transforms as a bidoublet of SO (4) ~ SU(2)l x SU(2)r and are identified as the SM Higgs 
doublet. The \x components of the remaining gauge fields have the following b.c: 



Wt ~ (+,+), B ~(+,+) 



(-,+), z> -(-,+), 



(A.2) 
(A.3) 



where a = 1, 2, 3, 6 = 1, 2, the L, indices correspond to the SU(2) l X 5 £7 (2)# decomposition 
of 50(4), and 



/J 



d + alx 



(A.4) 



Here we have defined X± = (X\±X2)/y/2, and X_ ^ satisfies (+, — ) b.c. Also, g$Xi = 9bX 2 = 
V2g 5X - 

The SM quarks are embedded in bulk fermions transforming in the fundamental repre- 
sentation of SO(5), 5 = (2,2) © (1, 1) with X + charge 2/3 and —1/3 for the up and down 



sectors, respectively. The odd fields couple as in Eq. (|2.3|). The b.c. are as follows: 



(2,2)£ 



(!,!)?(--) 



(2,2)| 



9ir( — ) 
(l,l)g(+,+) 



(2,2)«(+-) (2,2)£(-+) 
(l,l)t(-+) (!,!)&(+-) 



(2,2)i 2 



92L(++) 



(2,2 



,92 



<?2ii( ) 



(2,2)1 (+") (2,2)% (-+) 
(1,1)1 (-+) (1,1)^(+-) 



(A.5) 



;i,i)i 2 ( 



(M)S(+,+) 



(2,2)£(+-) (2,2)*(-+) 
(l,l)i(-+) (1, !)*(+-) 



We have displayed the field content according to their SO(4) decomposition, and + and 
— represent, respectively, Neumann or Dirichlet boundary conditions at the corresponding 
brane. Note that the choice of parities above seems to allow for two SM doublet zero modes 
per generation, coming from q\L and q^L- It is actually only the symmetric combination 
(QiL + Q2l)/V2 that has a zero mode, the odd combination being coupled to a UV-localized 
chiral fermion with a large mass. The 0(4) x U(l)x symmetry at the IR brane allows for 
the following mass mixing terms: 



m u (2,2^ 1 (2,2)^ + M tt (l,ir R 1 (l,l)^ + m d (2,2)^(2,2)^ + M,(l, 1)^(1, l)Uh.c. (A.6) 

Notice that the localized mass terms do not involve the odd fields in Thus, neglecting 
inter-generational mixing, we have a total of eight parameters per quark generation, four 



</2 , 
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bulk masses, denoted in units of the bulk curvature k by c qi , c q2 , c u and Cd (the exchange 
symmetry forces c u - = c u ), and four mass terms m U} d, M u d- The first two generations have 
a negligible effect on EWSB and, for the third generation, as long as c qi < 1/2, c qi < c q2 , the 
two multiplets related to the bottom (£g 2 ,£<f) do not play any significant role either. We are 
therefore left with four relevant parameters, c qi , c u , m u and M u . Note however that there 
is a phenomenological constraint on the value of q. The reason is that, due to the chosen 
boundary conditions, the component (1,1) becomes ultralight for < —1/2 [H], easily 
violating limits from direct production at the Tevatron |60| . 



B The Higgs Potential 



The one-loop Higgs potential in 5D warped models does not admit a simple analytic expres- 
sion, which is related to the difficulty in deriving an explicit form for the KK mass spectrum 
and Higgs interactions for the relevant fields. A useful tool to derive at least an implicit but 
relatively compact form for the potential is achieved by using the gauge-fixing outlined in 
|61j in an holographic approach [62]. In this way, the Higgs potential is simply obtained by 
a rotation of the 4D holographic fields only (see [61] for details). The Higgs potential in our 
model is actually the same as the one considered in [25] , the only difference being the presence 
of the Z2-odd SO(5) multiplet The latter field, interestingly enough, leads to a Higgs 
contribution which has opposite sign with respect to that of 5D fermions admitting chiral 
zero modes. Like gauge fields, they contribute positively to the Higgs potential, pushing the 
minimum towards zero. In order to simplify the discussion, let us consider a simple 5D model 
with SU(2) gauge symmetry broken to U(l) by boundary conditions. The contribution to 
the "Higgs" potential given by two 5D fermion doublets ^f± whose left-handed components 
satisfy the boundary conditions 



*1 



>!(—), 

modulo an irrelevant constant term, can be written as 

d 4 p 



V L (-+) 



V±(h) 



(27T)' 



log 



1 ±c 2/i n(p,m±) 



(B.l) 



(B.2) 



In Eq. (B.2), ± stands for the contributions of ^ + and ^~ , respectively, p is an Euclidean 
momentum and C2h = cos(2h/ fh). All the non-trivial information on the mass spectrum is 
encoded in the form factor H(p, m±), with p = ^Jp^p^ and m± the bulk 5D mass terms. 
If \H(p, m±)\ < 1 over the whole integration region, then one can expand the log term in 



Eq. (B.2), which clearly shows the opposite contribution to the potential given by the two 
5D fermion fields. One has 

\n -i 



n(p, to) 



G (c) _ G + _(c) 

G~&j ~ G^j 



G- 



+ 



Ga 



G_+(c) G ++ (c) 



(B.3) 
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where G± t ±(c) = G± t ±(c,p, z m , z vv ) are combination of the Bessel functions J and Y: 

G r , tV >(c,p,z 1 ,z 2 ) = J c + v i/2(iP Zi)Y c+v , 1/2 (ip z 2 ) - Y c+r/1/2 (ip z 1 )J c+ri , 1/2 (ip z 2 ) , (B.4) 

where rj,rj' = ±, and c = z vv m is the usual dimensionless mass term in warped space. It is 
not difficult to check that |II(p, m±)\ < 1, it exponentially vanishes when p — > oo and reaches 
its maximum value 1 just at the origin p = 0p^| 

Let us now turn to the actual Higgs potential in our model, which is complicated by 
the presence of localized IR mass terms and a non-trivial gauge symmetry breaking pattern. 
The most relevant fields contributing to the potential are the gauge fields and the SO(5) 
multiplets , f^, £ 9l and £ g2 . In the holographic approach, these contributions are encoded 
in the holographic gauge fields W and Z, the bottom and top quarks, and by a holographic 
Z2-odd fermion component of £~. 

The gauge field contribution reads (in the absence of any localized IR or UV terms) 



Vg{h) 



d* P 

WT 4 



2 log 1 + s 



. n_ - n_ 
' 211+ 



+ log 1 + 4- 



,n_ 



il 



2hi 



'l + cos26» 



w , 



where 9w is the SM weak mixing-angle and 



n± 



G T -(l/2) 
G T+ (l/2) 



(B.5) 



(B.6) 



in terms of the functions defined in Eq. (B.4). In Eq. (B.5), the first and second log terms 
correspond to the holographic W and Z contributions, respectively. 

Let us now consider the fermion contribution to the Higgs potential. Since by symmetry 
no IR mass terms can be introduced for £~, whose components satisfy boundary conditions of 
the (+—)/(—+) type, its contribution to the potential is given by V-(h) in Eq. (B.2), which 
can be rewritten (again, modulo constant terms) as 

d 4 p 



Vf(h) 



-2N r 



2tt 



log 



1 + 4 



2IIn 



(B.7) 



where N r 



3 is the QCD color factor and 

G-+(c u ) 



n, 



IT 



G — (c u ) ' 

G- + (Cn 

G 



_ + G++(c u ) 
(cu) G + -(c u ) 



(Bi 



The holographic bottom and top quark contributions are more involved. The top contribution 
can be written as 



Vto P (h) 



-2N r 



d 4 p 

Wf 4 



log 



2 n i 

1 + Sh zn% 



i + 4 



nf \ 

2Up 



22 (M?l_ 
h h 2U^U q Q 



(B.9) 



17 A similar result holds in flat space where Ii(p,m) — — (p 2 +m 2 )/(m 2 + p 2 cosh(2L\/p 2 + m 2 )), with L 
being the length of the segment. 
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in terms of the following form factors: 

N+ tU (l/M u ) 



TT C/ 



nf 

Mi 



D quu (l/M u ) 



N+ UU (1/M U ) N+ uU (m u ) 
D gitU (l/M u ) D qi 



+ 



D guu (m u ) D g2id (m d ) ' 
NzJmJ , N-(1/M u ) 



91 

D gitU (m u ) 

4 



+ 



11, 

D quu (l/M u ) ' 
m u l/M u 



D quU (m u ) D quU (l/M u ) 



(B.10) 



with 



7V±(M) = G- ± {c i )G+±{c j ) + M 2 G +± {c i )G- ± {c j ), 
Dij(M) = G-+{c i )G+-{c j ) + l^G+ + {c i )G--{c j ). 



(B.ll) 



In the above formulae, m u and M u are dimensionless IR brane mass terms (see appendix |A[) . 



The bottom quark contribution is also given by Eq. (B.9), provided one makes the substitu- 
tions c qi «-> c q2 , c u <-> Cd, m u <-> mrf, M u «-> in Eqs. (B.9) and (B.10). 

It is useful to consider various limits for the IR brane mass terms. When m u — > 0, 
M u — > 0, the form factor M" vanishes and the potential splits in the two contributions 



n n 
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111, with Ilo and III given in Eq. (B.8). In the opposite 



given by and £91,92- As expected, in this limit the contribution given by £;+ reduces to 
that of with 11^ 
limit, m u 

(++) / ( ) boundary conditions, and is given by Vjt(h) in Eq. (B.2 ). This is again expected, 



00, M u — ► 00, the contribution given by turns into that of a fermion with 



since the large IR brane mass terms effectively change the boundary conditions of the fermions 
on the IR brane. In the limits in which M u — * 00, m u — * or viceversa, IT" — > 0, M" — * 0, 
and hence does not contribute at all to the potential. As far as £91,92 are concerned, one 
should recall that bottom and top quarks are distributed in both multiplets £ 91 and £ g2 and 
that a localized UV right-handed fermion doublet is necessary to get rid of an unwanted zero 
mode left-handed doublet [25]. Such a localized field implies that the contributions of £ gi 
and £ 92 are entangled also in the limit of vanishing or very large IR brane terms. That is 



why IIq in Eq. (B.10) depends on both c qxm , c U) ^ and m u ^. Notice, however, that even in 
this case, if M u — ► 00, m u — > or viceversa, IT| X vanishes, so that the whole top contribution 
to the potential vanishes in this limit! Such result is particularly clear in the approach 
followed here, where the Higgs dependence of the holographic Lagrangian is obtained by 
a SO(5) rotation of the 4D holographic fields. The rotation gives rise to Higgs-dependent 
terms when SO(5) is broken at the IR brane by boundary conditions. In the limit M u — > 00, 
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m u — > or viceversa, the singlet and the bidoublet components of the SO(5) multiplets , 
£ qi and £ 92 have the same boundary conditions at the IR brane, implying SO(5) invariant 
boundary conditions and hence no couplings of such fermions with the Higgs. Actually, there 
is a whole one-dimensional family of boundary conditions, when m u = 1/M U , for which 



II" = n^ 1 = Mi = 0, as is clear from Eq. (B.10). Needless to say, similar results hold when 
considering the IR brane terms and M^. 

The total one-loop Higgs potential is finally given by 

V To t{h) = V g (h) + V top (h) + V bottom (h) + Vf(h) . (B.12) 
C The Trilinear Couplings 

The holographic approach, used in the previous appendix to derive the effective Higgs po- 
tential, can also be efficiently applied to the computation of the trilinear couplings c/l and 
gu of the DM candidate A"_ with the third generation quarks and the Z2-odd fermions. As 
explained in the main text, these couplings are an essential ingredient to compute the DM 
relic abundance. 

For simplicity, we neglect EWSB effects, which introduce only a slight change (< few %) 
in the computation of the DM relic abundance, an accuracy that we have no interest in 



achieving. Starting from Eqs. (3.4) and (3.5), our aim is to extract the trilinear couplings 



c/RtR fi-ijj- t R (for c u < 0) and g^L fi-ip-,L (for c u > 0) where ip—,R is the right-handed 
component of the NLOP singlet state and ip- t L is the left-handed component of the NLOP 
field in the bidoublet coupled to the top. In the holographic approach, g^/n are obtained by 
writing the bulk-to-boundary propagators for the holographic fields t L / R , ijj-^/R an d and 



integrating the vertices (3.4) and (3.5) over the compact space, taking care of normalization 
factors required to define canonically normalized fields. We can consider the cases c u > and 
c u < at the same time. The (tL,tn) system can be described by choosing as holographic 
degrees of freedom the UV values of the LH components of the £ ?1 (and £g 2 ) multiplet and 
the RH ones of the field (for m u ^ 0). In this way we obtain the quadratic holographic 
Lagrangian, 

c = t L i^t L+ t R t^ (c.i) 

P z vv P z uv 

with p = ^JPpP^ and a Minkowskian momentum and the form factors and Uq are as 



in Eq. (B.10) but with p — > —ip 



is 



The holographic fields are not mass eigenstates but rather a superposition of all the 
states in the 4D KK tower. As a consequence, we need some care to extract the couplings 
of the mass eigenstates from the holographic Lagrangian. To do this we use the fact that 
the holographic fields coincide with the (non-canonically normalized) mass eigenstates when 



3 The factor l/z m ensures that the holographic fields have canonical dimension in 4D. 
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we put them on mass-shell. This means that we can obtain the couplings of the KK states 
by simply computing the holographic couplings on-shell, provided we also find the correct 
normalization factors for the fields. 

For the top field, the wave functions of the relevant on-shell (i.e. p = 0) components along 
the multiplet are given by 




2 — r 
-17l u Z u Zt 

1 / Z \2+c 
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t L = f+(z)t L , 



(C.2) 



Z VV ^Z v y/ 



tR = fR~{z)tR- 



To find the correct normalization factors, Z^ R , we must require that the action have the 
canonical form for massless fields in a series expansion in p around the on-shell momentum. 
The L and R normalization factors are thus given by 



i n q ' u 



lim 

p^o z uv p 



i am> u 

hm 

p^o z VY op 



(C.3) 



Due to the absence of zero modes, the holographic description of the Z2-odd fields is more 
involved, since the action diverges for p — > 0. Nevertheless we can still use the holographic 
approach by a suitable expansion around the on-shell momentum p = m oc id, where m 0( id is a 
mass eigenvalue. 

The NLOP states that couple to and ti, ip-^ and are contained in the multiplet 

£~ as in Eq. (A. 5), namely £7^ = (1,1)^~ and the relevant component of (2,2)^", that we 
denote by £7. Both fields satisfy (-1 — ) b.c, so that we can choose as holographic fields the 
UV values of the corresponding components: tp R = yjz vv S,r( z vv) and ip^ = V z vv 
The holographic Lagrangian for this system at the quadratic level is 
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£ = Vr + Vl Vl ' 

p z vv p z vv 



(C.4) 



where Ul = — (IIo +IIi/2) 1 and 11^ = Ilo, evaluated again with Minkowski signature. The 
relevant on-shell field wave functions are 



1 



z \ 5 / 2 G — (c u ,-im R ,z m ,z) 
' G — (c u , —im R , z m , z V y) 



z vv ^ z vv 
1 / z \5/2 a 

Z VV ^ Z VV 



R 



++ {c u ,-tm L ,z m ,z) 
C++(c u , —im L , z m , z uv ) 



f R { Z )^R^ 

/r(*)V£- 



(C.5) 



19 In presence of EWSB, when the various field components are mixed by the Higgs VEV, it is better to 
choose holographic fields of the same chirality within a multiplet (see j6l]). In the present case the two 
components are independent of each other, so we can simply use holographic fields with different chiralities. 
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The normalization factors Z L R can be found by expanding the action around the on-shell 
momentum. One gets 



20 



{Z- LR f= lim lim _J_(^L 2 n ^ (C 6) 



where m LjR are the masses of the NLOP. 

The holographic description of the X_ gauge field is easily obtained in an axial-type gauge 
X_ ; 5 = 0, as described in {61]. The quadratic holographic Lagrangian for the transverse part 
of the X- gauge field is found to be 

c = -J—x t _ 4iP n.x t J t , (c.7) 

where X l _ = ^z vv X t i(z vv ) is the holographic field. The on-shell field wave function is 

Y t M _ * G ++ (l/2,-im x _,z m ,z) 9t _ 9t 

X -- {Z) ~ ^G ++ (l/2,^m x _,z 1R ,z vv ) X - = f*-W X -*> ( C - 8 ) 

and its normalization factor, , is 

Z J_ = lim -^^^ = lim AfH5=). (C.9) 
The interaction term can finally be written as 

= f<fcX?L/Jl r dz {^Y {fL/ R fX-fi/ R ) + h - C -> ( C - 10 ) 

where gsx is the 5-D £/(l)x gauge coupling. Taking into account the normalization factors, 
we finally get: 

From the above relations it is possible to get an upper bound on gi and gn by noticing that 
both couplings are maximal when tn and ti reside mostly in the £ u multiplet, which occurs 
for M u — > oo and m u — > oo, respectively]^] When M u — > oo and c u ~ — 1/2, or m u — ► oo and 

20 The factor 1/2 which appears in these expressions is due to the fact that the fermions we are describing are 
massive. In the holographic approach the two chiralities of the fermions (say 4>l,r) are described by the same 

holographic field (4>l), so that on-shell 4>l = 4>t and (f>R = ]6/p4>L- In this case cj>(j/—m)(f> — > (f) L 2(p — rn)^)/p<pL- 
21 Notice that in the limit M u — > oo, m u — > oo, the components of the f u multiplet related to it and 
in acquire (++) b.c, while the corresponding £ qi components satisfy (H — ) boundary conditions [see also 
paragraph after Eq. (B.lllj. 
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c u ~ 1/2 one gets: 

/ \ 3/2 

{Z~ R f ~ (Z L ) 2 ~ (Z X _) 2 , ~ / L (z) * — fx.(z) , 



1 



(Z+) 2 ~ bg(W^Dv) . /£(*) - -7^= f — ) . (C12) 

/, \l-2csn / \V2-Cm / \3/2 

(^ + ) 2 -^( — ) l0g(W%v), 



3/2 



^uv V ^uv / V ^uv 



Using Eqs. (C.12), we see that the integral over the wave functions in Eq. (C.ll) reduces to 



the normalization condition for the gauge field both for gn and §l, giving simply 

,- ,2 95X 2 

IflL/fll < o / i ; , N = o9x , (C.13) 
o V 2; uvlog(z IR ./2 U v) 3 

where in the latter equality we have used the relation between the 5D coupling g§x and the 



4D coupling gx- This proves Eq. (3.6) in the main text 
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